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ABSTRACT 

The analysis of laminar mixed convection in boundary layers 
adjacent to a vertical, continuously stretching sheet in the presence of 
radiation and heat sink has been studied. It has been assumed that thermal 
conductivity is varying linearly with temperature. The stretching velocity, 
surface temperature and surface concentration varying linearly with the 
distance from stagnation point. The governing flow equations have been 
transformed into set of coupled ordinary differential equations using 
similarity transformation. These equations have been solved by using 
Runge-Kutta method with shooting technique. The effects of Prandtl 
number, Schmidt number, sheet stretching parameter, heat source or sink 
parameter and thermal conductivity parameter, on the boundary layer 
velocity, temperature and concentration profiles have been discussed in 
detail for both the cases of assisting and opposing flows. The results 
show that magnitude of temperature increases with increasing values of 
thermal conductivity parameter and velocity and temperature both 
increase with increasing values of heat sink parameter. Also magnitude 
of heat transfer coefficient decreases with increasing thermal conductivity 
and with heat sink parameter for assisting as well as opposing flow. 

Key words : Stretching sheet, variable thermal conductivity, 
heat source, radiation, stagnation point, shooting technique. 



INTRODUCTION important problem in engineering processes 

with applications in industries such as in hot rolling, 
The flow over a stretching sheet is an wire drawing, crystal growing, continuous casting, 
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glass fiber production, extrusion of plastics and 
paper production. The study of heat and mass 
transfer is necessary for determining the quality 
of the final products of such processes. 

The flow near a stagnation point has 
attracted many investigators during the past several 
years. In 1988 Karve and Jaluria presented 
flow and mixed convection transport from a 
moving plate in rolling and extrusion processes 
and in 1991 they done the numerical simulation 
of thermal transport with a sheet in material 
processing. Mahapatra and Gupta studied the 
heat transfer in the steady two-dimensional 
stagnation point flow of an incompressible fluid 
over a stretched sheet with a velocity proportional 
to the distance from the stagnation point, also 
with effect of viscous dissipation of the fluid in 
2002. Mixed convection boundary layer flow 
in the stagnation point flow toward a stretching 
sheet is presented by Ishak in 2006. The boundary 
layer flow for an incompressible viscous fluid 
near a stagnation point at a heated stretching sheet 
placed in a porous medium, using Lie-group method 
is analyzed by Boutros. In 2007 El-Aziz studied 
the heat and mass transfer of electrically 
conducting fluid having temperature dependent 
viscosity and thermal conductivity along a 
stretching sheet with the effect of Ohmic heating. 
In 2007 Mahapatra and Dholey studied the 
magneto hydrodynamic fluid flow towards a 
stretching sheet. In 2009 Jing ZHU found the 
analytical solution to a stagnation point flow and 
heat transfer over a stretching sheet using 
Homotopy analysis. Recently in 2009 D. Pal 
studied effect of thermal radiation on heat and 
mass transfer stagnation point flow towards a 
stretching sheet. 

In all the previous studies in mixed 
convection flow towards a stretching sheet, the 



effect of heat source with radiation and variable 
thermal conductivity is not studied. Aim of this 
paper is to study the effects of variable thermal 
conductivity and heat source on steady two 
dimensional mixed convection stagnation point 
flow towards a stretching sheet in the presence 
of radiation effect. It is assumed that the stretching 
velocity, surface temperature and surface 
concentration vary linearly with the distance 
from stagnation point. The governing fluid flow 
equations are transformed into non dimensional 
form using similarity transformation and solved 
by using Runge- Kutta method with shooting 
technique. 

Mathematical Formulation of Problem : 

The mixed convection, two-dimensional, 
laminar, boundary layer flow due to a stretching 
vertical heated sheet in a viscous, incompressible 
and electrically conducting fluid in the presence 
of heat source and radiation with variable thermal 
conductivity has been made. The stationary 
coordinate system has its origin located at the 
centre of the sheet with the x-axis extending 
along the sheet, while the y-axis is measured normal 
to the surface of the sheet and in the positive 
direction from the sheet to the fluid. The coordinate 
system has been shown in adjacent figure. The 
stretching velocity u„(x), surface temperature 
T w (x) and surface concentration C w (x) vary 
linearly with the distance from stagnation point 
O. The governing fluid flow equations for the 
problem under consideration are: 

Equation of continuity, 
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Figure (a) 



Equation of momentum, 



fluid, Tqo an d Cqq are temperature and concen- 
tration at outer side of boundary layer, k' is 
variable thermal conductivity of fluid, Q is heat 
source parameter and D is mass diffusion 
coefficient. In equation (2) "+" sign is corres- 
ponding to assisting flow and "-"sign refers to 
an opposing flow. The stretching surface has 
temperature T w , concentration C„ and the free 
stream temperature and concentration are T^ 
and Cqo respectively, with T„> T^ and C„> 
Cqo- In this study the radiation heat flux in the 
x-direction has also been considered. Using 
Rosseland diffusion approximation the radiation 
heat flux is given by: 



4a* dT 4 

3k* dy 



(5) 



Keeping the temperature difference within the 
flow is assumed to be sufficiently small, so that 
T 4 may be expressed as a linear function of 
temperature T, i.e. 



^ + »l^ + v^ ±gW T-T») ±8 MC-C») 
ox ay ox dy 2 

(2) 

Equation of heat transfer, 



T 4 =4T3T-3T4 



dx dy pC p dy 



dyh pC p dy pC p 



Equation of mass transfer, 



dC dC n d 2 C 
u — + v — = D 

dx dy dz 2 



(3) 



(4) 



where u and v are velocity component along x 
and y directions, g is acceleration due to gravity, 
T and C are temperature and concentration of 
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The initial and boundary conditions are: 

u= u„(x) =U 0 x, v=0, T=T„(x) =Too+ bx, 
C=C„(x) =Coo+ d x, at y=0 

u=u e (x)=ax, T-^Tqo, C^Cqq, as y^co (7) 

Here Uo, a, b, d are constants. u e is velocity of 
the flow external to the boundary layer. 
The similarity transformations used are: 



(8) 



V = (U 0 v)^xf(T 1 ),T 1 = [l^| y 
9v|/ 



d\i> 
u = — ,v 
dy 



dx 



and k'=k(l+e9) 
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T-T C-C 

9(T1)=^— ^-,4>(T1)=> 



T - T 



c - c 



here \|/ is stream function, r| is similarity 
variable, 8 is thermal conductivity parameter, 
k is thermal conductivity of fluid, 0 and § are 
dimensionless temperature and concentration. 
Using above transformations equations (2) - 
(4) are converted into following set of ordinary 
differential equations: 



(9) 



f m -HM-l -f 2 +ff"±},e±o(t) = o 



Sc(ff-f '(|>) + ())" = o 



(10) 

(ii) 



where Sc is Schmidt number, Ec is Eckert 
number, X is thermal buoyancy parameter, 8 is 
solutal buoyancy parameter. 

The transformed boundary conditions are: 

f'(0)=l, f(0)=0, 0(0)=land <|)(0)=1 at r|=0 
f=-jf-, 0=0 and (|)=0 as r|-»oo (12) 



Where « = —,Re x = ^and 9 W = ^,Rd = ^» 
pC p v k k* 



G r c = gPc(Cw-C.)x3 = gP T (T ; -T.)x3 

are local Grashof number and local solutal 
Grashof number respectively. Also A. = ^ and 



o = are thermal and solutal buoyancy 
parameters. The physical quantities of practical 



interest are skin friction coefficient Cf, the local 
Nusselt number Nu x and the local Sherwood 
number Sh x . These quantities are given as 



xq„ 



C f =— ,Nu x 

pu 2 w k'(T w -T„) D(C W -C„) 



,Sh x 



(13) 



Here 



q« =- 



lit 



-f 



y=0 

H6a*T 3 



3k : 



+ k' 



<9T 



\ dy 



y=0 
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(1+80)0" + f Rd [1+(0 W -1)0] 3 0" + e0' 2 + 
4Rd[l+(0 w -l)0] 2 0 ,2 (0 w -l)+Prtf0'-f0'+S0] =0 1* 



Now after using (9) and (14), we obtain 



C f Re x /2 =f"(0), 

Sh x 



Nu x 



/Re x 



Re x 



(l+e+fRde 3 w )e'(0) 
(T+s8) 



(15) 



RESULTS AND DISUCUSSION 



The Runge-Kutta method with shooting 
technique has been used to solve the set of coupled 
ordinary differential equations Ar|=0.01 and 
1100=8. These results show good agreement 
with Pal (2008) for a/U 0 =l, Pr=0.72, 0„=1.1, 
X=l, Rd=l, 5=0.5, Sc=0.5, S=0 and 8=0. The 
value of skin friction coefficient f '(0) in the 
paper by Pal is 0.64268 for assisting flow and 
in the present study the value is of f '(0) for 
same set of parameters is 0.642959. The value 
of dimensionless velocity f '(r|), dimensionless 
temperature 0(r|) and dimensionless concen- 
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tration ®(r|) have been computed for different 
values of dimensionless parameters and shown 
from figures (1) - (10). These results are observed 
for both assisting as well as for opposing flow. 
Figure (1) presents the temperature profile with 
varying values of thermal conductivity parameter 
(s) for assisting and opposing flow. It has been 
observed that magnitude of temperature 
increases with increasing values of thermal 
conductivity parameter. Figures (2) and (3) 
present the effect of Prandtl number on velocity 
and temperature of fluid. These figures show that 
magnitude of velocity and temperature both 
decreases with increasing values of Prandtl 
number (Pr=0.7, 5, 10). In figures (4) and (5) 
effect of heat source parameter on velocity and 
temperature have been shown. It has been seen 
that magnitude of velocity and temperature 



both increase with increasing values of heat 
source parameter. The effects of Schmidt 
number on velocity and concentration profiles 
have been shown in figures (6) and (7). These 
figures show that magnitude of velocity as well 
as the magnitude of concentration both 
decrease with increasing values of Schmidt 
number (Sc=0.5, 2, 5). From figures (8) - (10) 
effect of stretching velocity parameter (a/Uo 
=0.5, 1, 1.5, 2) on velocity, temperature and 
concentration of fluid have been presented. In 
figure (8) it has been seen that magnitude of 
velocity increases with increasing values of 
stretching velocity parameter. In figures (9) and 
(10), it has been seen that magnitude of tempe- 
rature and concentration both decrease with 
increasing values of stretching velocity para- 
meter. 
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Figure 1. Temperature profile with thermal 

conductivity parameter 
(Pr=7, a/U 0 =l, S=0.1, Rd=0.5, X=-l, 5=0.5, 
Sc=0.5, 0 W =1) 



Figure 2. Velocity profile with Prandtl 
number 

(8=1, a/U 0 =l, S=0.1, Rd=0.5, X=-l, 5=0.5, 

sc=o.5, e„=i) 
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Figure 3. Temperature profile with prandtl 
number 

(8=1, a/U 0 =l, S=0.1, Rd=0.5, X=-l, 5=0.5, 

sc=o.5, e w =i) 



Figure 4. Velocity profile with heat source parameter 

(8=1, a/U 0 =l, Pr=7, Rd=0.5, X=-l, 5=0.5, 

sc=o.5, e w =i) 
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Figure 5. Temperature profile with heat 

source parameter 
(s=l, a/U 0 =l, Pr=7, Rd=0.5, A=-l, 5=0.5, 

sc=o.5, e„=i) 



Figure 6. Velocity profile with Schmidt 
number 

(s=l, a/U 0 =l, Pr=7, Rd=0.5, X=-l, 5=0.5, 

s=o.i,e„=i) 
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Figure 7. Concentration profile with Schmidt 
number 

(s=l, a/U 0 =l, Pr=7, Rd=0.5, X=-l, 5=0.5, 

s=o.i,e w =i) 



Figure 8. Velocity profile with Strerching 

velocity parameter 
(8=1, Sc=0.5, Pr=7, Rd=0.5, X=-l, 6=0.5, 

s=o.i,e„=i) 
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Figure 9. Temperature profile with 
Strerching velocity parameter 
(s=l, Sc=0.5, Pr=7, Rd=0.5, X=-l, 8=0.5, 

s=o.i,e„=i) 



Figure 10. Concentration profile with 
Strerching velocity parameter 
(s=l, Sc=0.5, Pr=7, Rd=0.5, X=-l, 5=0.5, 

s=o.i,e„=i) 
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The values of skin friction coefficient, heat transfer coefficient and mass transfer coefficient 
for assisting and opposing flow have been given in table 1 and 2, for different values of 
dimensionless parameters. 



Table 1. Values of skin friction coefficient, heat transfer coefficient and mass transfer 
coefficient with different values of dimensionless parameters for assisting flow 



X=l, 5=0.5, e w =l.l,Rd=0.5 


Pr 


Sc 


S 


a/U ( , 


s 


f'(0) 


9'(0) 


f(0) 


0.7 


0.5 


0.1 






0.650039 


-0.58543 


-0.93701 


7 


1 


0.1 






0.476752 


-1.51852 


-0.91494 


7 


5 


0.1 






0.41772 


-1.78382 


-1.28459 


7 


0.5 


0.1 






0.352132 


-1.77387 


-2.84539 


7 


f 0.5 


0.5 






0.463501 


-1.55798 


-0.91347 


7 


0.5 


1 






0.486871 


-1.22499 


-0.91543 


7 


0.5 


0.1 


0.5 




-0.161044 


-1.68675 


-0.81306 


7 


0.5 


0.1 


2 




2.389276 


-2.01820 


-1.09576 


7 


0.5 


0.1 


1 


2 


0.470360 


-1.49734 


-0.91407 



Table 2. Values of skin friction coefficient, heat transfer coefficient and mass transfer 
coefficient with different values of dimensionless parameters for opposing flow 



1=1, 5=0.5, 9 w =l.l,Rd=0.5 


Pr 


Sc 


S 


a/Uo 


e 


f'(0) 


9'(0) 


f(0) 


7 


0.5 


0.1 






-0.476331 


-1.68752 


-0.85585 


7 


5.0 


0.1 






-0.366038 


-1.70709 


-2.75523 


0.7 


0.5 


0.1 






-0.73360 


-0.50304 


-0.81615 


10 


0.5 


0.1 






-0.443503 


-2.02828 


-0.85906 


7 


0.5 


0 






-0.471969 


-1.74698 


-0.85625 


7 


0.5 


1.0 






-0.533382 


-1.02939 


-0.85023 


7 


0.5 


0.1 






1.636983 


-1.96508 


-1.06844 


7 


0.5 


0.1 




0 


-0.444990 


-2.22646 


-0.85848 


7 


0.5 


0.1 




2 


-0.501742 


-1.40201 


-0.85339 
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In this study, it has been seen that 
magnitude of skin friction coefficient increases 
with thermal conductivity parameter, heat source 
parameter and also with stretching velocity 
parameter for both assisting and opposing flow. 
But magnitude of skin friction coefficient decreases 
with increasing values of Prandtl number and 
Schmidt number for both observed cases. Further 
it has been noticed that magnitude of heat 
transfer coefficient decreases with increasing 
thermal conductivity parameter and with heat 
source parameter for assisting as well as 
opposing flow. It has also been observed that 
with increasing Prandtl number, Schmidt 
number and with stretching velocity parameter, 
the magnitude of heat transfer coefficient 
increases. In this study magnitude of mass 
transfer coefficient increases with stretching 
velocity parameter and with Schmidt number for 
both observed cases and there is no significant 
change in heat transfer coefficient with other 
parameters. 
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